Diffusional growth of wetting droplets 
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The diffusional growth of wetting droplets on the boundary wall of a semi-infinite system is 
considered in different regions of a first-order wetting phase diagram. In a quasistationary ap- 
proximation of the concentration field, a general growth equation is established on the basis of 
a generalized Gibbs-Thomson relation which includes the van der Waals interaction between the 
droplet and the wall. Asymptotic scaling solutions of these equations are found in the partial-, 
complete- and pre-wetting regimes. 
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The physics of wetting phenomena has attracted much interest in recent years, both from experimental and 
from theoretical [§]-[)) points of view. Whereas initially static properties dominated the discussion, the interest has 
shifted more recently to the dynamics of wetting P,[l0|--|l7|| . In many experimental situations the formation of a wetting 
layer starts with the nucleation of droplets on the boundary wall of the system. The central question therefore is the 
temporal evolution of the droplet profile. 

There are essentially two different types of dynamic behavior of a liquid surface droplet. The first is a spreading 
process which e.g. dominates, if a droplet of a non- volatile liquid is overheated from below to above a wetting transition 
point. Such processes are driven by hydrodynamic modes of the liquid, and they have extensively been discussed in 

■ the literature ]qjlfj|-|l^ . The second mechanism is the phase transformation (condensation or evaporation) between 
the liquid and the vapor phase of the droplet. This is driven by particle diffusion in the vapor, and e.g. is the 

f — ' dominating process in the growth of supercritical droplets in a metastable situation. 

Whereas the diffusional growth of a homogeneous wetting layer has been discussed in the literature , this seems 
not to be the case for surface droplets. The present paper deals with the diffusional growth of a supercritical droplet 
from a super-saturated vapor. This process is accompanied by the creation of latent heat, and it will be assumed that 
heat transport as well as other hydrodynamic modes are fast compared to the diffusion. As a consequence the droplet 

OA . is isothermal and always has a shape which minimizes its free energy at a given volume. The time dependence of this 

■ shape is the main object of interest in this paper. 

The excess free energy of a wetting film of local thickness f(x) on a planar boundary wall of a semi- infinite system 
H ■ can be written in the form ||l8|,[l9| 



n h [f] = I d 2 x l(vff + v(f)-hf , (i) 
o 

where 7 is the interface stiffness, h is the difference of the chemical potential from that of the saturated vapor and 
, V(f) is an effective interface potential. The field h can be expressed by the difference between the vapor concentration 
c and its value c at saturation, so that in linear order 



c = c (1 + Th) . (2) 

The form of the potential V(/) corresponding to a first-order wetting transition is sketched in Fig. 0. There, for 
temperatures T less than the wetting temperature T w , the global minimum of V(f) is at / = f , whereas for T > T w 
this minimum becomes metastable in favor of the global minimum at diverging film thickness. For / — » 00 we assume 
V oc f x ~~ a , where a — 3 for nonretarded and a — 4 for retarded van der Waals interactions ||. 

Homogeneous (i.e. fix) = const) minima of the excess free energy TLh[f\ (i-e. the global minima of V{f) — hf) 
determine the phase diagram, shown in Fig. ||. In the region h > 0, where the liquid bulk phase is stable, a film 
of infinite thickness forms on the wall in thermal equilibrium. On the line h = 0, which means bulk coexistence of 
the liquid and vapor phases, the first-order wetting transition occurs at T = T w , where f — fo for h = — 0, T < T w 
(partial wetting), and f — 00 for T > T w (complete wetting). From the transition point a prewetting line hpiT) 
extends into the region h < where the vapor phase is stable in the bulk. This line separates a region (below h p (T)) 
where the wall is covered by a thin film from a region (/i > hpiT)) where the wall is covered by a thick film. The 
jump in film thickness along the prewetting line vanishes at the prewetting critical point T pw . The partial wetting 
line h = 0, T < T w and the prewetting line h p iT) together form a first order line concerning the wetting properties 
of the system Q . 

If the system is quenched from below to above the first order line (for example by increasing the pressure), the 
phase transition is initialized by the formation of critical droplets on the wall (provided, one stays within the surface 
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spinodal lines, shown in Fig. ||). The shape of these droplets is qualitatively different in different regions of the phase 
diagram [E2hM|. Axisymmetric profiles f(r) can be calculated via the saddle point equation 



5H h /6f(r) = (3) 
with the natural boundary conditions of a droplet profile 

/'(0) = 0, lim f(r) = f . (4) 

r — >oo 

As illustrated in Fig. EL this leads to spherical (in the squared gradient approximation of Eq. (jl]) parabolic) caps in 
the partial wetting regime, to flat cylindrical droplets (pancakes) in the prewetting regime l| and to ellipsoid-like 
droplets in the complete wetting regime |2^]. The saddle point SHh / Sf(r) — has an unstable growth mode but the 
volume preserving shape fluctuations are stable ■ 

Assuming that the volume growth of the droplet is slow, the diffusion in the surrounding concentration field c 
becomes quasistationary and can be approximated by the Laplace equation 

D Ac = 0, (5) 

where D is the diffusion constant. At far distances from the droplet the concentration field is given by the system 
concentration Coo(i) which is time-dependent in a supersaturated system (h > 0) because of the phase-separation 
process in the metastable bulk phase. The normal derivative of the concentration field on the boundary wall of the 
system vanishes because there is no diffusion flux into the wall, i.e. the Neumann boundary condition 

D d x c\^ = (6) 

has to be fulfilled. 

To obtain a well defined diffusion problem the boundary condition on the surface as well as the actual shape of the 
supercritical droplet need to be specified. Motivated by the slow diffusional growth of the droplet the concentration 
field close to the droplet surface is assumed to be in local thermal equilibrium. Therefore the local chemical potential 
h(x) at the droplet surface is given by h(x) = SHo / Sf(x) which due to (||) corresponds to a concentration 

*(*) = co(i + r^), (7) 

which can be denoted as a generalized Gibbs-Thomson relation for wetting droplets. The expression STCq / Sf consists 
of a term 7 times the local curvature K of the droplet interface plus an interaction term dV / df. Neglect of the 
interaction term reduces (0) to the classical Gibbs-Thomson relation c s — cq(1+AK) with the capillary length A = Tj. 
It identifies the concentration at a curved interface as the concentration cq for a flat interface modified by a linear 
curvature correction. 

The assumption of fast hydrodynamic modes (compared with the diffusional growth) implies that the shape of 
a growing droplet can be calculated by minimizing its free energy under the constraint of a fixed droplet volume 
£l(t). Technically, this variational calculation leads again to Eq. (g) with the boundary conditions (^) but now with 
h in (|J) replaced by a function /in(t) which includes a Lagrange parameter corresponding to O(t). Consequently 
the growing supercritical droplet always looks like a critical droplet at a different time-dependent chemical potential. 
With increasing volume Q(t) the corresponding field /i^m approaches the first order line, where eventually the volume 
of the droplet diverges. In this sense wetting droplets grow along isotherms towards the first order line, as illustrated 
in Fig. g. 

The saddle point equation (|3j) with the fixed volume constraint is equivalent to STto / Sf = /lo(t)- 

Via Eq. @ this 

implies that the Dirichlet boundary condition of the constrained equilibrium droplet is given by 

c s (t) =c (l + Th m ) (8) 

and therefore independent of x. Especially in the complete wetting or prewetting case, where the droplets are not 
spherical, one would expect a non trivial boundary condition having the classical Gibbs-Thomson condition in mind. 
Additionally corrections due to the potential V(f), which determine the shape of the droplets in these regions, have 
to be taken into account. Nevertheless both effects add up in such a way that Eq. (0) can be written as Eq. (||) for a 
droplet in a volume-constraint equilibrium showing that c s is constant along the droplets surface! 

Now, the ^-independent Dirichlet boundary condition (^) allows to use an electrostatic analogy to solve the qua- 
sistationary diffusion problem (^)-(R) for the growing droplet ^8j. To fulfill the Neumann condition (||), the system 
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(including the droplet) is mirrored at the boundary wall of the system. Then the field AnDc is identified with an 
electric potential which also obeys the Laplace equation. The normal derivative of the field, i.e. the diffusion flux 
density on the droplet surface field, corresponds to the charge density of a conductor with the shape of the droplet 
including its mirror image. Consequently, the total volume growth of the droplet corresponds to the total charge, 
which is given by the capacity C of the conductor times the potential difference between the surface and infinity. This 
ultimately leads to the droplet growth equation 

tl = 4<7rDC(t) [ Coo (t) - c s (i)] = 4irDTC(t) [h(t) - h m ] (9) 

where C depends on the droplets profile and therefore is implicitly time-dependent The difference [h(t) — /ifi(t)] 
may be interpreted as the supersaturation of the system with respect to the droplet. 

Eq. (||) together with Eqs. ([!]), (||) and (||) allow to determine selfconsistently the growing droplet profile if h(t) is 
known. For a given volume f2, the droplet profile can be calculated from by Eqs. (||) and (^) with a conveniently 
chosen Lagrange multiplier hn . Then the capacity C of the conductor represented by the droplet plus its mirror image 
is calculated. Insertion of C, Iiq and the chemical potential h into Eq. @ yields the droplet growth rate and 
integration of (^|) eventually determines £l(t). 

In practice, for large droplets, the calculation can be facilitated by the use of scaling properties of critical droplets 



close to the first-order transition line 1 22 |23J . At temperatures T > T w the wetting droplets on the wall nucleate either 
as ellipsoid-like droplets at h > (complete wetting) or as pancake like droplets at h p (T) < h < (prewetting). In 
both cases the droplets grow along an isotherm towards the prewetting line (/in(t-»oa) ~> hp(T))- This means that 
they eventually become pancake-like droplets with a constant height but diverging radius R(t) so that the capacity 
of large droplets is given by the capacity of a flat disk C(t) oc R(t). In the case where the initial quench leads to a 
supersaturated bulk system (h > 0) the volume will phase separate until it reaches h = whereas for initial values 
h < the vapor bulk phase is stable and h remains constant in time. In either situation the difference [h — /in] 
approaches a non-vanishing constant so that Eq. (^) yields tl oc SI 1 / 2 or il cx t 2 which implies 

R<xt, (10) 

and only is determined by the time dependent capacity, i.e. the increasing diffusive coupling to the environment. 

Wetting droplets at T = T w in a supersaturated system (h > 0) also are not spherical. Their radius R scales as 
R cx /iq~ < - <t+1 ' ) / 2<T , their central height F as F cx /in 1/,<T , and consequently their volume as £1 cx hn~ <,a+2 ^ a [fj2| . 
Therefore, the profile of a growing wetting droplet becomes flatter and approaches a disk with capacity C cx R cx 
hn~( a+1 ' 1 2a . In a supersaturated system (i.e. h > 0) there are not only wetting droplets on the wall, but also 
droplets in the bulk. The set of growing bulk droplets reduces the supersaturation in a Lifshitz-Slyozov- Wagner type 
way as h cx i -1 / 3 M-E7|. With this input the wetting droplet growth equation (0) can be written as 



2(g+l) _£±i . I/O 

h n /in - oc /in 2(7 At l ' i 



(11) 



for large droplets. This leads to the asymptotic growth law hn oc t 3 ("+ 3 ) which due to the above scaling properties 
for R and F implies 

2(g + l) 4 

R oc t^+3j , p oc 1 3( CT+3 ) . (12) 

Finally, in the partial wetting regime T < T w , h > 0, the wetting droplets are spherical caps, and therefore their 
growth properties are similar to those of bulk droplets, i.e. 

R(xt 1/3 (13) 

The evaluation of the difference [h — /in] in Eq. (JoJ) can only be done in a theory where diffusional interactions 
between surface and bulk droplets are taken into account pi| . At late stages partial wetting droplets in systems with 
a temperature corresponding to a contact angle > tt/2 will shrink because [h — ha] turns negative for each droplet, 
whereas droplets at a temperature with < ir/2 will grow. 

One of the basic ingredients of the present calculation is the Neumann boundary condition (^) . It derives from the 
fact that there is no diffusion flux through the wall. Even if the wall locally is in a non-wet state it always is covered 
by a film of microscopic thickness /q. If somewhere in such a region Eq. (^) were not valid, the film would thicken 
there and the interface would run out of the microscopic minimum of the interface potential shown in Fig. [l]. Due 
to the generalized Gibbs-Thomson relation (^), which includes a term dV / df, the local concentration on top of the 
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surface would then increase and the film would evaporize until it again reaches the former height fo- Thus, up to 
fluctuations, Eq. (^) will be valid. 

The calculation of the supercritical droplet shape is based on the assumption of fast hydrodynamic modes compared 
to the droplets diffusional growth. This assumption may become questionable, if a prewetting droplet becomes very 
large. However at this very late stage the coalescence of different droplets will be dominant anyway. 
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FIG. 1. Sketch of an effective interface potential that shows a first order wetting transition as T is raised from T < T w to 
T > T w . 
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FIG. 2. Wetting phase diagram: The first order line of wetting transitions consisting of the partial wetting line h — 0, T < T w 
and the prewetting line h = h p (T) is marked by a solid line. The different types of growing wetting droplets along isotherms are 
shown in the different regions of the phase diagram. The dotted lines refer to the surface spinodals which enclose the nucleation 
regime. 
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